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We separate the contributions from the secondary star and 
from other sources with the aim of determining a photometric 
distance. Simple model calculations show that the accretion 
stream and the cyclotron source on the accreting white dwarf 
are minor contributors to the infrared light, allowing an accu- 
rate determination of spectral type and absolute flux of the sec- 
ondary star. The unilluminated backside of the secondary star 
as seen in eclipse has spectral type dMO to dM0+. Its dered- 
dened magnitude is K — 13.58 at orbital phase = (eclipse). 
Using the calibrated surface brightness of M-stars and the pub- 
lished mass of the secondary, M2 = 0.46 M , we obtain a dis- 

1 /2 

tance d = 600 + 25 pc which scales as MJ . The radius of the 
Roche-lobe filling secondary exceeds the main-sequence radius 
of an MO star by 21 +1 f j%. The debated origin of the infrared 
^? light of V1309Ori has been settled in favor of the secondary 
star as the main contributor and an accurate distance has been 
fvq [ derived that will place estimates of the luminosity and synchro- 
1 ■ nization time scale on a more secure basis. 
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Abstract. The first phase-resolved JHK light curves of the eclipsing polar (AM Herculis binary) VI 309 Ori are presented and 
interpreted. 



1. Introduction 

The soft X-ray source RXJ05 15.6+0105 was discovered in 
the ROSAT All-Sky-Survey (Beuermann & Thomas 1993), 
quickly identified as an eclipsing polar (or AM Herculis bi- 
nary) (Garnavich et al. 1994; Shafter et al. 1995; Walter et 
al. 1995), and named VI 309 Ori. It turned out to be a key 
object for studies of the synchronization of polars because of 
its long orbital period of 7.98 h, twice as long as that of any 
other polar, and the implied larger separation of the compo- 
nents. At the same time, the magnetic field strength of about 
60 MG (Garnavich et al. 1994; Shafter et al. 1995) and the ac- 
cretion rate of ( 1 .0 + 0.5) X 1 0~ 9 M yr~ 1 (Beuermann & Burwitz 
1995;Harrop-Allinetal. 1997; Schwarz et al. 2005) are within 
the ranges typical of long-period polars. Luminosity-based es- 
timates of the accretion rate depend sensitively on the bolo- 
metric correction, the mass of the primary, and on the assumed 
distance. Published distances of VI 309 Ori differ considerably, 
ranging from about 500 pc (Garnavich et al. 1994; Staude et 
al. 2001) to greater than 1500pc (Harrop-Allin et al. 1997). 
The synchronization of a polar requires that the magnetostatic 
torque between the assumed magnetic dipole moments of the 
two stars in the system overcomes the accretion torque on the 
white dwarf (King et al. 1994; Wu et al. 1994), with the lat- 
ter depending directly on the accretion rate. The parameters 
derived so far place V 1 309 Ori near the limit of being synchro- 
nized and measuring its luminosity and mass transfer rate is 
important for understanding the physics of polars in general. 

In this paper, we derive an accurate distance from the 
brightness of the secondary in VI 309 Ori. This result is 
based on our previous calibration of the surface brightness 



* Based on observations obtained under run ID 01BK6 at the 
Canada-France-Hawaii Telescope (CFHT) which is operated by the 
National Research Council of Canada, the Institut National des 
Sciences de l'Univers of the Centre National de la Recherche 
Scientifique of France, and the University of Hawaii. 



of M-dwarfs in standard and non-standard photometric bands 
(Beuermann 2000) and will place estimates of the luminosity, 
accretion rate, and synchronization time scale on a more secure 
basis. 

2. Observations and Data Analysis 

VI 309 Ori was observed with the CFHT-IR camera at the 
CFHT 3.6-m telescope on Mauna Kea, Hawaii, for 6.7 hours on 
26/27 October 2001 . During the observations, the star was in its 
normal bright state. The filters employed were J, H, K', and K. 
Almost complete phase coverage was obtained with a gap only 
near orbital phase cf> = 0.4, where <p — indicates the center 
of the eclipse for the ephemeris of Staude et al. (2001). Five 
15 sec exposures taken in dithering mode in a single filter were 
averaged to yield an individual data point. Exposures in J, H, 
and K were obtained consecutively and the sequence repeated 
in turn, leading to a time resolution of 13 min. Variability on 
shorter time scales is lost and only in this sense are the observa- 
tions in the three bands quasi-simultaneous and not truly simul- 
taneous. UKIRT faint standards FS 150, 152, 153, and SA92- 
342 were observed as standard stars and the results reduced into 
the standard Mauna Kea infrared system (Tokunaga & Vacca 
2005). We interprete these data in combination with our earlier 
BVRI photometry and low resolution optical spectrophotome- 
try (Shafter et al. 1995) and the photometry of Katajainen et 
al. (2003). 

3. Light curves 

Figure 1 shows the light curves of VI 309 Ori in the bands 
J, H, and K folded over the orbital period of 7.98 h using 
the ephemeris of Staude et al. (2001). They display a quasi- 
sinusoidal modulation with the eclipse superimposed. The 
depth of the eclipse is 0.31 mag, 0.20mag, and 0.15 mag in 
J, H, and K, respectively, much less than in the optical where 
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it reaches 2.1 mag in B and nearly 4 mag in U (Shafter et al. 
1995; Katajainen et al. 2003). The light outside eclipse arises 
from the accretion stream and and the cyclotron source in the 
accretion spot on the white dwarf. This combined source evi- 
dently still contributes in the infrared, although much less than 
in the optical. Given our low time resolution, we can not mea- 
sure ingress and egress times of the eclipsed source. Our data, 
however, are consistent with the eclipse being total for a time 
interval similar to that observed in the optical (~ 40 min). 

At first glance, the double-humped structure of the light 
curve seems to suggest that we see the ellipsoidal light modula- 
tion of the secondary star. However, that is refuted by the phas- 
ing of the primary eclipse and the finite eclipse depth. Many 
polars (and non-magnetic CVs) contain light sources which 
mimic the ellipsoidal modulation. In the case of polars, this 
can be the accretion stream which contributes maximum flux 
when seen sideways and reaches a minimum when the plane of 
the accreting field line crosses the line of sight. In VI 309 Ori, 
this is the case at ^ 0.91. A double-humped light source 
can also originate from cyclotron beaming and the varying flux 
level as the angle between the accreting field line at the white 
dwarf and the line of sight varies. We find that the maxima at 
</> = 0.2 and <p = 0.7 are of nearly equal height in the infrared, 
while the <p = 0.2 maximum is the dominant one in the optical. 
The available results combined suggest that the light curves are 
shaped by the varying geometrical aspect, mutual occultation 
of the stream by itself or the stream and the secondary star, and 
wavelength-dependent illumination by the soft X-ray and cy- 
clotron sources on the white dwarf. We refrain from modeling 
the light curves in detail. 

We wish to derive the spectral energy distributions (SED) 
inside and outside eclipse and, for this purpose, we have to 
model the ellipsoidal light modulation of the secondary star. 
We assume that the eclipse is total (Schmidt & Stockman 2001 ; 
Staude et al. 2001) and that the flux observed at eclipse center 
arises from the backside hemisphere of the secondary, i.e. the 
one opposite to the Lj point. We calculate the ellipsoidal light 
curves using Nightfall 1 . This program uses the Roche geome- 
try, the limb darkening coefficients of Claret (1998), and the 
gravity darkening coefficients of Claret (2000) that describe 
the transition from high-temperature radiative stars to low- 
temperature convective stars. It allows the use of Hauschildt's 
model atmospheres for non-irradiated stars, but for a simple ap- 
plication including irradiation one fares better with the black- 
body approximation. We have adjusted the resulting ellipsoidal 
light curves at </> — and show them both, for the case of 
vanishing irradiation (solid curves) and for irradiation by a 
source at the location of the white dwarf with a luminosity of 
10 34 erg cm _2 s _1 (dashed curves). This luminosity corresponds 
to the rather high accretion rate of 1 .5 x \0 9 M yr~' . The influ- 
ence of irradiation is seen to be moderate in the infrared, while 
the minimum at <p — 0.5 nearly vanishes in the optical for the 
parameters chosen here. We obtain the SEDs at </> = and at 
(p = 0.21 by combining our infrared photometry with the our 
earlier BVRI photometry, our 3600-9 140 A spectrophotometry 
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Fig. 1. Quasi-simultaneous near-infrared light curves of VI 309 Ori in 
the J, H, and K bands folded with the binary orbital period. Also shown 
are the calculated light curves of the ellipsoidal modulation of the sec- 
ondary added. The K-band light curve has been shifted upwards by 
0.5 mag in order to avoid overlap. The solid (dashed) curves refer to 
the unirradiated (irradiated) secondary star. The data sets at the bot- 
tom of the figure show the difference, representing infrared light from 
sources other than the secondary (see text). 

(Shafter et al. 1995), and the photometry of Katajainen et al. 
(2003). 

4. Spectral energy distributions 

Fi gure 2 shows the SED at cp — 0, dereddened with £b-v — 
0.15 (Schmidt & Stockman 2001). The Shafter et al. spec- 
trophotometry in eclipse is slightly adjusted by a quadratic 
function in wavelength to fit the BVRI photometry from 
the same reference. The data point for the U-band is from 
Katajainen et al. (2003). The SED is characterized by a color 
V - K — 3.64 (see Tab. 1). This value closely resembles that of 
a dwarf of spectral type M0 to M0+. The M0 SED is shown as 
a dashed line in Fig. 2 that is adjusted to the observed VIJHK 
data points and is seen to fit extremely well. The SED of an Ml 
dwarf (dotted line) is clearly too red. In the U-band, an excess 
of the photometry over the M0 SED may exist and possibly also 

Table 1. Dereddened magnitudes of the secondary star in VI 309 Ori at 
orbital phase tp = (center of eclipse) compared with the magnitudes 
of M0 and Ml field adjusted to the same mean flux level as VI 309 Ori. 



1 http://www.hs.uni-hamburg.de/DE/Ins/Per/Wichmann/Nightfall.html 
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Fig. 2. The spectral energy distribution of the remnant light in the cen- 
tral phase intervals of the eclipse is consistent with an origin from the 
secondary in V1309Ori alone. Besides our JHK photometry (filled 
circles), published UBVRI photometry (open circles) and the Shafter 
et al. (1995) spectrophotometry (solid line) are shown. Included for 
comparison is the flux distribution of a young-disk M0 dwarf (dashed 
line) adjusted to the mean flux level of V1309 Ori and that of an Ml 
dwarf adjusted at V (dotted line). 



Fig. 3. Spectral energy distribution of the light from sources other than 
the secondary at orbital phase <p = 0.21. The symbols are as in Fig. 2. 
The (subtracted) contribution of the secondary at <f> = 0.21 is shown as 
the dashed line. A spectrum representative of the hydrogen emission 
of the stream (solid line) and a cyclotron spectrum from the accretion 
spot (dot-dashed line) are shown for comparison. Their normalization 
is arbitrary (see text). 



a slight excess in B. A possible cause is an atmospheric struc- 
ture of the secondary star in VI 309 Ori which differs slightly 
from that of a field star, even on the non-irradiated hemisphere. 
Alternatively, a 2% leak from the bright phase spectrum would 
also explain the UB excess, but the weak emission lines ex- 
pected in that case are not visible in eclipse. Furthermore, 
the Hubble Space Telescope phase-resolved ultraviolet spec- 
troscopy (Schmidt & Stockman 2001; Staude et al. 2001) and 
the model calculations of Katajainen et al. (2003) suggest that 
the eclipse of the accretion stream is total. No deviation from 
the SED of an M0 to M0+ star beyond the ~ 0.03 mag scat- 
ter is seen for all bands from V to K. The excellent agreement 
confirms that the infrared flux observed in eclipse represents 
the secondary star. 

The SED for <p = 0.21, with the contribution from the sec- 
ondary star (dashed line labeled 'Sec') subtracted, is shown 
in Fig. 3. The dereddened Shafter et al. spectrophotometry is 
shown but not adjusted to the photometry in this case. Although 
the optical and infrared observations are not simultaneous, an 
internally consistent distribution emerges which falls off with 
wavelength approximately as fx °c AT 26 . The secondary star is 
obviously the dominant source in the infrared, while the other 
sources dominate in UBVR and both contribute equally in the 
I band. The secondary contributes 3/4 of the observed flux in 
K which clearly contradicts the claim by Harrop-Allin et al. 
(1997) of a contribution of only 20%. These authors based their 
statement on the lack of detectable spectroscopic M-star fea- 
tures in their infrared spectrum. Although we have no spectral 
resolution in the infrared, we consider our conclusion on the 
SED of the secondary star firm and suggest instead that the 
Harrop-Allin et al. infrared spectrum, which was taken around 



(j) = 0.5, is affected by irradiation of the secondary star and may 
differ from that of a field star. 

Although we do not embark on a complete modeling of the 
<p = 0.21 SED, we comment semi-quantitatively on the two 
main contributing light sources. The Shafter et al. spectropho- 
tometry at this phase shows the Balmer and the Paschen jumps 
in emission along with the appropriate emission lines suggest- 
ing a temperature around 1 0000-20000 K. Our heuristic model 
consists of two main contributions to the observed SED. As 
solid curve labeled 'Str', we show the emission from an isother- 
mal hydrogen plasma at a temperature of 15000K with two 
subcomponents of different density and geometrical extent: the 
first subcomponent has a particle density of 5 x 10 13 ctrT 3 , a 
depth of 10 9 cm 2 , unity relative projected area, and represents 
the outer accretion stream; the second subcomponent has a par- 
ticle density of 3 x 10 14 crrT 3 , a depth of 3.5 x 10 8 crrT 2 , a rel- 
ative projected area of 0.35, and represents the inner accretion 
stream a factor of two in radius closer to the white dwarf. Stark 
broadening is included and the resolution of the resulting spec- 
trum adjusted to the 35A FWHM resolution of the spectropho- 
tometry. Only the sum of the two subcomponents is shown with 
an arbitrary normalization. This spectrum fits the shape of the 
observed optical spectrum reasonably well and is seen to fall off 
in the infrared with a similar slope as the observed SED. As the 
dot-dashed curve labeled 'Cyc', we show the cyclotron emis- 
sion from an isothermal plasma with a temperature of 5keV 
in a magnetic field of 63 MG for a combination of two depth 
parameters, log A = 6 which yields a largely optically thick 
spectrum, and log A = 2 which accounts for the optically thin 
emission lines. The latter fits the two faint and broad bumps at 
about 5800A and 8500A reported by Shafter et al. (1995) that 
represent the second and third cyclotron harmonics in a field 
of 63 MG. Our model spectrum demonstrates that the fourth 
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harmonic near 4300A is already weak in the log A = 2 spec- 
trum, while the first harmonic at 1 .7/im is largely optically thick 
(and off the scale). Hence, we do not expect to observe more 
than the two detected cyclotron harmonics. The high density 
of the accretion stream at the white dwarf deposits much of 
the accretion energy at sub-photospheric levels (Walter et al. 
1995; Schwarz et al. 2005) and explains why cyclotron line 
emission and polarization are weak (Garnavich et al. 1994; 
Shafter et al. 1995; Buckley & Shafter 1995). The observed 
overall SED is quite well reproduced by a linear combination 
fx = 2.5 f a + 0.5 / cyc . The white dwarf represents a minor 
contribution as shown by the thick line labeled 'WD' which 
is based on the upper limit to the white dwarf contribution as 
derived by Staude et al. (2001). A quantitative model could 
be constructed along these lines, but the optical-depth effects 
are substantial and one may expect a sizeable fraction of the 
stream emission to fall in the EUV (Lyman jump in emission), 
where no observations exist. It is doubtful, therefore, whether 
convincing results for the bolometric luminosity and the accre- 
tion rate can be obtained this way; we do not attempt it here. 

5. Distance of V1 309 Ori 

We determine the distance to VI 309 Ori by the surface bright- 
ness method (Barnes & Evans 1976; Bailey 1981; Ramseyer 
1994; Beuermann et al. 1999; Beuermann 2000) which con- 
nects the distance d, the stellar radius R, the surface brightness 
Fx at wavelength A, and the dereddened observed flux fx by the 
inverse square law 

This relation is often quoted with the surface brightness and the 
flux expressed in magnitudes (S x and nix) as 

log d = (mx - S x)/5 + 1 + log(R/R Q ). (2) 

Using a sample of main sequence stars with appropriate met- 
alicity, the F A or S x for different photometric bands have been 
calibrated as functions of color or spectral type (Beuermann et 
al. 1999; Beuermann 2000, 2006). The distance is found to be 
proportional to the stellar radius of the secondary star which is 
given by Roche geometry. 

Staude et al. (2001) quote the mass range of the sec- 
ondary in VI 309 Ori as 0.4-0.6 M , with a most probable value 
of M2 = 0.46 M G . The radius of the Roche-lobe filling sec- 
ondary scales as M^ 3 and is nearly independent of the white 
dwarf mass. The appropriate projected area is that of the Roche 
lobe as seen at <p = at an inclination of 78° (Staude et 
al. 2001). The corresponding effective mean radius is R = 
4.88 x 10 10 cm = 0.701 R . We use Eq. (2) with K = 13.58 
and S K = 3.91 for V - K = 3.64 (Beuermann 2000) to obtain 
d = 602 pc. If there is some blue excess in the SED of the sec- 
ondary it might be more appropriate to use V - K larger by 0. 10 
and Sk = 3.94 and which yields d = 594 pc. The statistical and 
systematic errors in the measurement of K and the calibration 
of Sk of about 0.05 each lead to a 4% error in d. We settle on 
d = 600 ± 25 ^7 P c > where the second error refers to the full 



range of M2 from Staude et al. (2001). Instead of K and 5k, we 
might have used any other photometric band from V to H. The 
results would be similar because the secondary corresponds so 
closely to the SED of an M0 to M0+ field star. The error, how- 
ever, is smallest in K because Sk(V - K) has the largest lever 
arm and the shallowest slope of all surface brightness vs. color 
relations. As an additional example, we quote the result for the 
narrow-band surface brightness F7500 at 7500A which was also 
calibrated by Beuermann (2000). From the SED in Fig. 2 one 
finds a dereddened flux /7500 = 7.19 x 10" 16 ergcm _2 s _1 A _1 . 
With F7500 = (1.16 + 0.07) x lO^rgcnrV'A- 1 (Beuermann 
2000) for a star of spectral type M0 to M0+, Eq. (1) yields 
d = 635 pc with a slightly larger error than above 2 . 

Our result is in excellent agreement with the distance ob- 
tained from the comparison of the brightness of the secondary 
star in VI 309 Ori with that of the mean component of the bi- 
nary YY Gem which yields 625 pc (Staude et al. 2001, cor- 
rected to the Hipparcos distance of YY Gem by Schwarz et al. 
2005). The much larger distance suggested by Harrop-Allin et 
al. (1997), d > 1500pc, based on the assignment of 80% of 
the K-bwd flux to the accretion stream, can be excluded. 

The distance modulus to V1309Ori is K - M K = 8.89+° 
and with a 4n average magnitude of the unheated secondary 
of K = 13.48, its absolute TT-band magnitude would be Mr = 
4.59+Q -\g, where the errors refer to the full range of secondary 
masses of 0.4 to 0.6 M , with 0.46 M the nominal mass. A 
dM0+ main sequence star of near-solar metalicity has Mk = 
5.0 and a radius R given by log(/?/R ) = -0.22 (Beuermann et 
al. 1998, 1999), while the volume radius of the Roche-lobe fill- 
ing secondary in V1309Ori is R = (5.04+°g) x 10 10 cm with 
log(fl/R ) = -0.14+°;°* Hence, we find that the Roche-lobe fill- 
ing secondary is expanded over a main sequence star of the 
same spectral type by (21 +1 i)%, as noted first by (Garnavich 
et al. 1994). This oversize may indicate moderate nuclear evo- 
lution. On the other hand, the secondary in V471 Tau which 
does not fill its Roche lobe has a radius 18% larger than a 
main sequence K-star, probably because a substantial fraction 
of its surface is covered by star spots (O'Brien et al. 2001). 
Activity may, therefore, be an alternative cause for the ob- 
served inflation of the secondary in VI 309 Ori. A measurement 
of the metal abundances supplemented by evolutionary model- 
ing may possibly distinguish between the two possibilities. 

6. Conclusions 

We have presented the first phase-resolved infrared photome- 
try of VI 309 Ori and deduced the spectral energy distributions 
(SED) of the secondary star and of the additional light sources 
outside eclipse for the 3600-22000 A range. The secondary is 
of spectral type M0 to M0+ and contributes 75% of the light 
in the /f-band. We derive a distance of 600 pc from its /f-band 
magnitude. The SED at orbital phase (f> = 0.21 can be plausibly 
explained as the superposition of emission from the accretion 
stream and the cyclotron source on the white dwarf. 

2 The distance of 745 pc quoted by Beuermann (2000) did not in- 
clude dereddening and was based on an a slightly different calibration 

Of F 7500- 



6 



K. Reinsch et al.: The distance of V1309 Ori 



Acknowledgements. This research was supported by the Korea 
Astronomy Observatory and Space Science Institute Research Fund 
2002 and by the Deutsches Zentrum fur Luft- und Raumfahrt e.V. 
(DLR) under grant 50 OR 99 03 6. We are deeply indebted to Y.- 
B. Jeon for generous support of our observing proposal within the 
KAO-CFHT collaboration. We thank Boris T. Gansicke for providing 
his program to calculate the emission from a hydrogen plasma and 
Frederic V. Hessman for comments on the manuscript. 

References 

Bailey, J. 1981, MNRAS, 197, 31 
Barnes T.G., Evans D.S. 1976, MNRAS, 174, 489 
Beuermann, K. 2000, New Astr. Rev., 44, 93 
Beuermann, K. 2006, A&A, submitted 

Beuermann K., Baraffe I., Kolb, U., Weichhold, M. 1998, 
A&A, 339,518 

Beuermann K., Baraffe I., Hauschildt, P. 1999, A&A, 348, 524 
Beuermann, K., & Burwitz, V. 1995, ASP Conf. Ser., 85, 99 
Beuermann, K., & Thomas, H. C. 1993, Adv. Space Res. 
13(12), 115 

Buckley, D. A. H., & Shafter, A. W. 1995, MNRAS, 275, L61 
O'Brien, M. S., Bond, H. E., & Sion, E. M. 2001, ApJ, 563, 
971 

Claret, A. 1998, A&AS, 131, 395 
Claret, A. 2000, A&A, 359, 289 

Garnavich, P. M., Szkody, P., Robb, R. M., Zurek, D. R., & 

Hoard, D. W. 1994, ApJ, 435, L141 
Harrop-Allin, M. K., Cropper, M., Potter, S. B., Dhillon, V. S., 

& Howell, S. B. 1997, MNRAS, 288, 1033 
Katajainen, S., Piirola, V., Ramsay, G., et al. 2003, MNRAS, 

340, 1 

King, A. R., Kolb, U., de Kool, M., & Ritter, H. 1994, MNRAS, 
269, 907 

Ramseyer, T. F. 1994, ApJ, 425, 243 
Schmidt, G. D., & Stockman, H. S. 2001, ApJ, 548, 410 
Schwarz, R., Reinsch, K., Beuermann, K., & Burwitz, V. 2005, 
A&A, 442, 271 

Shafter, A. W., Reinsch, K., Beuermann, K. et al. 1995, ApJ, 
443,319 

Staude, A., Schwope, A. D., & Schwarz, R. 2005, A&A, 374, 
588 

Tokunaga, A. T., & Vacca, W. D. 2005, PASP, 1 17, 421 
Walter, F. M., Wolk, S. J., & Adams, N. R. 1995, ApJ, 440, 834 
Wu, K., Wickramasinghe, D. T., Bailey, J., Tennant, A. F. 1994, 
PAS A, 11, 198 



